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Abstract. Arnold, Falk, and Winther recently showed [Bull. Arner. Math. 
Soc. 47 (2010), 281-354] that linear, mixed variational problems, and their 
numerical approximation by mixed finite element methods, can be studied 
using the powerful, abstract language of Hilbert complexes. In another recent 
article [arXiv: 1005.4455], we extended the Arnold-Falk-Winther framework by 
analyzing variational crimes (a la Strang) on Hilbert complexes. In particular, 
this gave a treatment of finite clement exterior calculus on manifolds, general- 
izing techniques from surface finite element methods and recovering earlier a 
priori estimates for the Laplace-Beltrami operator on 2- and 3-surfaces, due 
to Dziuk [Lecture Notes in Math., vol. 1357 (1988), 142-155] and later Demlow 
[SIAM J. Numer. Anal, 47 (2009), 805-827], as special cases. In the present 
article, we extend the Hilbert complex framework in a second distinct direction: 
to the study of semilinear mixed problems. We do this, first, by introducing 
an operator-theoretic reformulation of the linear mixed problem, so that the 
semilinear problem can be expressed as an abstract Hammcrstein equation. 
This allows us to obtain, for semilinear problems, a priori solution estimates 
and error estimates that reduce to the Arnold— Falk— Winther results in the 
linear case. We also consider the impact of variational crimes, extending the 
results of our previous article to these semilinear problems. As an immediate 
application, this new framework allows for mixed finite element methods to be 
applied to semilinear problems on surfaces. 
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1. Introduction 

The goal of this paper is to extend the abstract Hilbert complex framework of 
Arnold, Falk, and Winther [5j — which they introduced to analyze certain linear 
mixed variational problems and their numerical approximation by mixed finite 
elements — to a class of semilinear mixed variational problems. Additionally, we 
aim to analyze variational crimes in this semilinear setting, extending our earlier 
analysis of the linear case in Hoist and Stern [24 . 

1.1. Background. Briining and Lesch [5] originally studied Hilbert complexes as a 
way to generalize certain properties of elliptic complexes, particularly the Hodge 
decomposition and other aspects of Hodge theory. More recently, Arnold, Falk, and 
Winther [4] showed that Hilbert complexes are also a convenient abstract setting 
for mixed variational problems and their numerical approximation by mixed finite 
element methods, providing the foundation of a framework called finite element 
exterior calculus (see also [3]). This line of research is the culmination of several 
decades of work on mixed finite element methods, which have long been used with 
great success in computational electromagnetics, and which were more recently 
discovered to have surprising connections with the calculus of exterior differential 
forms, including de Rham cohomology and Hodge theory [SJ H7J [551 121] • For this 
reason, Hilbert complexes are a natural fit for abstract methods of this type. 
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Another recent development in this area has been the analysis of "variational 
crimes" on Hilbert complexes (Hoist and Stern [23] )• By analogy with Strang's 
lemmas for variational crimes on Hilbert spaces, this work extended the estimates 
of Arnold, Falk, and Winther [3] to problems where certain conditions on the 
discretization have been violated. This framework also allowed for a generalization 
of several results in the field of surface finite element methods, where a curved 
domain is not triangulated exactly, but is only approximated by, e.g., piecewise 
linear or isoparametric elements. This research area was initiated with the 1988 
article of Dziuk |T7] (see also Nedelec (25]), with growing activity in the 1990s [T8"1[T2"] 
and a substantial expansion beginning around 2001 [2*21 HTl H§1 IH | |2"0" 1 H5 1 H5 1 H5] . 

Our main motivation for extending the estimates of Arnold, Falk, and Winther 
[3] and of Hoist and Stern [55], from linear to semilinear problems, is to enable 
the use of finite element exterior calculus for nonlinear problems on hypersurfaces, 
allowing for a complete analysis of the additional errors due to nonlinearity, as well 
as those due to surface approximation. 

1.2. Organization of the paper. The remainder of the article is structured as 
follows. In |Section 2| we give a quick overview of abstract Hilbert complexes and 
their properties, before introducing the Hodge Laplacian and the linear mixed 
problem associated with it. We then discuss the numerical approximation of 
solutions to this problem, summarizing some of the key results of Arnold, Falk, and 
Winther [4] on approximation by subcomplexes, and those of Hoist and Stern [24] 
on variational crimes. In |Section 3[ we introduce an alternative, operator-theoretic 
formalism for the linear problem, which — while equivalent to the mixed variational 
formulation — allows for a more natural extension to semilinear problems, due to its 
monotonia ty properties. We then introduce a class of semilinear problems — which 
can be expressed in the form of certain nonlinear operator equations, called abstract 
Hammerstein equations — prove the well-poscdness of these problems, and establish 
solution estimates under various assumptions on the nonlinear part. In |Section 4[ 
we extend the a priori error estimates of Arnold, Falk, and Winther 4 from linear 



problems to the semilinear problems introduced in Section 3 including improved 
estimates subject to additional compactness and continuity assumptions. Finally, we 
generalize the linear variational crimes framework of |24] to this class of semilinear 
problems. These last results allow the linear a priori estimates, established in [53] 
for surface finite element methods using differential forms on hypersurfaces, to be 
extended to semilinear problems. 

2. Review of Hilbert complexes and linear mixed problems 

We begin, in this section, by quickly recalling the basic objects of interest — Hilbert 
complexes and the abstract Hodge Laplacian — along with the solution theory for 
linear mixed problems in this setting. This provides the background and preparation 
for semilinear problems, which will be discussed in the subsequent sections. The 
treatment of this background material will be necessarily brief; we will primarily 
follow the approach of Arnold, Falk, and Winther [4], to which the interested reader 
should refer for more detail^ At the end of the section, we will also summarize 



1 This is largely a condensed version of the background material given in Hoist and Stern 1241 
Section 2] , from which we quote freely. We include it here in the interest of keeping the present 
paper self-contained, since the semilinear theory will depend, to a large degree, on several properties 
and results that have recently been established for the linear problem. 
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the results from Hoist and Stern [24j . analyzing variational crimes for the linear 
problem, in preparation for extending these results to the semilinear case. 

2.1. Basic definitions. First, let us introduce the objects of study, Hilbert com- 
plexes, and their morphisms. 

Definition 2.1. A Hilbert complex (W, d) consists of a sequence of Hilbert spaces 
W k , along with closed, densely-defined linear maps d fc : V k C W k — > V k+1 C W k+1 , 
possibly unbounded, such that d k o d fe_1 = for each k. 

■ ■ ■ > v k ~ x > v k > v k+1 > ■ ■ ■ 

This Hilbert complex is said to be bounded if d k is a bounded linear map from W k 
to W k+1 for each k, i.e., (W, d) is a cochain complex in the category of Hilbert 
spaces. It is said to be closed if the image d k V k is closed in W k+1 for each k. 

Definition 2.2. Given two Hilbert complexes, (W, d) and (W-^d'), a morphism 
of Hilbert complexes f:W—t W consists of a sequence of bounded linear maps 
f k : W k -> W' k such that f k V k C V' k and d' k f k = f k+1 d k for each k. That is, the 
following diagram commutes: 

■ ■ ■ > v k > V k+1 > ■ ■ ■ 

f k f k+1 

d' k 

■ ■ ■ > v' k > y>k+i > ■ ■ ■ 

By analogy with cochain complexes, it is possible to define notions of cocycles, 
coboundaries, and harmonic forms for Hilbert complexes. (This also gives rise to a 
cohomology theory for Hilbert complexes.) 

Definition 2.3. Given a Hilbert complex (W, d), the space of k-cocycles is the 
kernel 3 fe = kerd fe , the space of k- coboundaries is the image = d k ~ 1 V k ~ 1 , and 
the kth harmonic space is the intersection $) k = 3 fe H < B k± . 

In general, the differentials d fe of a Hilbert complex may be unbounded linear 
maps. However, given an arbitrary Hilbert complex (W^d), it is always possible to 
construct a bounded complex having the same domains and maps, as follows. 

Definition 2.4. Given a Hilbert complex (W, d), the domain complex (V, d) consists 
of the domains V k C W k , endowed with the graph inner product 

(u,v}y k = (u,v) wk + (d k u,d k v) wk+1 . 

Remark 1. Since d fe is a closed map, each V is closed with respect to the norm 
induced by the graph inner product. Also, each map d k is bounded, since 

\\d k v\\ vk+1 = \\d k v\\ wk+1 < \\v\\ wk + \\d k v\\ wk+1 = \\v\\ vk . 

Thus, the domain complex is a bounded Hilbert complex; moreover, it is a closed 
complex if and only if (W, d) is closed. 

Example 2.5. Perhaps the most important example of a Hilbert complex arises 
from the de Rham complex (f2(M),d) of smooth differential forms on an oriented, 
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compact, Riemannian manifold A/, where d is the exterior derivative. Given two 
smooth fc-forms u,v £ il fe (M), the L 2 -inner product is defined by 

( u ^ v )l^q(m) = uA*v= Mi 
Jm jm 

where * is the Hodge star operator associated to the Riemannian metric, ((•, •)) is the 
metric itself, and \x is the Riemannian volume form. The Hilbert space L 2 £l k (M) is 
then defined, for each k, to be the completion of Vl k (M) with respect to the L 2 -inner 
product. One can also define weak exterior derivatives d k : HQ, (M) C L 2 Q k (M) — ► 
HQ k+1 {M) C L 2 Q k+1 (M); the domain complex (HQ(M),d), with the graph inner 
product 

( u > w /Hn(M) = ( u > v )L 2 n{M) + (d«jdv)£2fi(M) , 

is analogous to a Sobolev space of differential forms. (For example, in M 3 , the 
domain complex corresponds to the spaces H 1 , H (curl), and H (div).) Finally, we 
mention the fact that both the L 2 - and iT-de Rham complexes are closed. For a 
detailed treatment of these complexes, and their many applications, see Arnold, 
Falk, and Winther [4]. 

For the remainder of the paper, we will follow the simplified notation used by 
Arnold, Falk, and Winther 0]: the W -inner product and norm will be written simply 
as (•,•) and ||-||, without subscripts, while the V- inner product and norm will be 
written explicitly as (v)y an d IHly 

2.2. Hodge decomposition and the Poincare inequality. For L 2 differential 
forms, the Hodge decomposition states that any fc-form can be written as a direct 
sum of exact, coexact, and harmonic components. (In M 3 , this corresponds to the 
Helmholtz decomposition of vector fields.) In fact, this can be generalized to give a 
Hodge decomposition for arbitrary Hilbert complexes; this immediately gives rise 
to an abstract version of the Poincare inequality, which is crucial to much of the 
analysis in Arnold, Falk, and Winther 0]. 

Following Bruning and Lesch [S], we can decompose each space W k in terms of 
orthogonal subspaces, 

w k = 3 fe e 3 k±w = 3 k n (b* © <s k1 -) © 3 fe±w = B* © fi k © 5 k±w , 

where the final expression is known as the weak Hodge decomposition. For the 
domain complex (V,d), the spaces 3 fe , and Sj k are the same as for (W, d), and 
consequently we get the decomposition 

where 3 k±v = ^ k±lv n V k . In particular, if (W, d) is a closed Hilbert complex, then 
the image 25 fc is a closed subspace, so we have the strong Hodge decomposition 

W k = <B k ffl-f^ ®3 k±w , 

and likewise for the domain complex, 

V k = ( B fc ffijo fc ffi3' £ - Lv . 

From here on, following the notation of Arnold, Falk, and Winther [J], we will 
simply write 3 fcJ ~ in place of 3 k±v when there can be no confusion. 
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Lemma 2.6 (abstract Poincare inequality). If (V,d) is a bounded, closed Hilbert 
complex, then there exists a constant cp such that 

\\v\\ v <c P \\d k v\\ v , V«e3 fcl 

Proof. The map d k is a bounded bijection from 3 fc± to *B' £+1 , which are both closed 
subspaces, so the result follows immediately by applying Banach's bounded inverse 
theorem. □ 

Corollary 2.7. If (V,d) is the domain complex of a closed Hilbert complex (W,d), 
then 

IMIv <cp||d fe v ||, Vi;e3 fe± . 

We close this subsection by defining the dual complex of a Hilbert complex, and 
recalling how the Hodge decomposition can be interpreted in terms of this complex. 

Definition 2.8. Given a Hilbert complex (W, d), the dual complex (W* , d*) consists 
of the spaces W k = W k , and adjoint operators d* k = (d fc_1 )* : V k * C W k -> F fe *_ x C 

■ ■ ■ < v k u «_^_ v k r ^ v k * +1 < 

Remark 2. Since the arrows in the dual complex point in the opposite direction, 
this is a Hilbert chain complex rather than a cochain complex. (The chain property 
d* k o d* k+1 = follows immediately from the cochain property d fe o d*^ 1 = 0.) 
Accordingly, we can define the k-cycles 3fc = kerd£ = *B feJ - w and k-boundaries 
Q3|! = d* k+l V k *. The fcth harmonic space can then be rewritten as S) k = 3 fc H 3£; 
we also have 3 fe = w , and thus 3 feJ_w = 25£. Therefore, the weak Hodge 
decomposition can be written as 

w k = W®^ k ®^i, 

and in particular, for a closed Hilbert complex, the strong Hodge decomposition 
now becomes 

W k = s B fc ©io fc ©<B^. 

2.3. The abstract Hodge Laplacian and mixed variational problem. The 

abstract Hodge Laplacian is the operator L = dd* + d*d, which is an unbounded 
operator W k — > W k with domain 

D L = {ue V k n V k * | du g V k * +1 , d*u G V 1 *- 1 } . 

This is a generalization of the Hodge Laplacian for differential forms, which itself is 
a generalization of the usual scalar and vector Laplacian operators on domains in 
K n (as well as of the Laplace-Beltrami operator on Riemannian manifolds). 
If u G Dl solves Lu = /, then it satisfies the variational principle 

(du,dv) + (d*u,d*v) = (f,v) , Vw G V k nV k . 

However, as noted by Arnold, Falk, and Winther [3], there are some difficulties in 
using this variational principle for a finite element approximation. First, it may be 
difficult to construct finite elements for the space V k <~)V k * . A second concern is the 
well-posedness of the problem. If we take any harmonic test function v G f) k , then 
the left-hand side vanishes, so (/, v) = 0; hence, a solution only exists if / _L $rj k . 
Furthermore, for any q G fj k — 3 fc H 3fc, we have dq = and d*q — 0; therefore, if u 
is a solution, then so is u + q. 
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To avoid these existence and uniqueness issues, one instead defines the following 
mixed variational problem: Find (a, u,p) € V k ~ 1 x V k x Sj satisfying 

(ff,T> - (u,dr) =0, VreF fe - 1 , 

(1) (do-, u) + (du, d«) + (p, v) = (/, «) , Vu € y fe , 

(u,?)=0, Vggij*. 

Here, the first equation implies that a = d*u, which weakly enforces the condition 
u € V k H Vj(?. Next, the second equation incorporates the additional term (p,v), 
which allows for solutions to exist even when / JL Sj k . Finally, the third equation fixes 
the issue of non-uniqueness by requiring u _L S) k . The following result establishes 
the well-posedness of the problem ([lj . 

Theorem 2.9 (Arnold, Falk, and Winther g], Theorem 3.1). Let (W,d) be a closed 
Hilbert complex with domain complex ( V, d) . The mixed formulation of the abstract 
Hodge Laplacian is well-posed. That is, for any f g W k , there exists a unique 
(a,u,p) g V k ~ x x V k x Sj k satisfying Moreover, 

\W\\v + \\u\\ v + \\p\\<c\\f\\, 

where c is a constant depending only on the Poincare constant cp in \Lemma 2.6\ 

To prove this, they observe that |T]) can be rewritten as a standard variational 
problem — i.e., one having the form B (x, y) — F(y) — on the space V k " 1 x V k x S) k , 
by defining the bilinear form 

B (a, u, p;r,v,q) = (a, r) — (u, dr) + (da, v) + (du, dv) + (p, v) — (u, q) 

and the functional F (r, V, q) = (/, v). The well-posedness of the mixed problem then 
follows by establishing the inf-sup condition for the bilinear form B (•, •) gl Theorem 
3.2], which shows that it defines a linear homeomorphism. This well-posedness 
result implies the existence of a bounded solution operator K : W k — > W k defined 
by Kf = u. 

2.4. Approximation by a subcomplex. In order to obtain approximate numeri- 
cal solutions to the mixed variational problem ([lj, Arnold, Falk, and Winther [4] 
suppose that one is given a (finite-dimensional) subcomplex Vh C V of the domain 
complex: that is, Vt C V k is a Hilbert subspace for each k, and the inclusion 
mapping ih : Vh c -s- V is a morphism of Hilbert complexes. By analogy with the 
Galerkin method, one can then consider the mixed variational problem on the 
subcomplex: Find {crh,Uh,Ph) € x V, k x Sj k satisfying 

(a h , T )-(u h ,dT)=0, VrgV^" 1 , 

(2) (da h ,v) + (du h , dv) + (p h ,v) = (/, v), Vve V k , 

(u h ,q)=0, VqeSj k h . 

For the error analysis of this method, one more crucial assumption must be 
made: that there exists some Hilbert complex "projection" tt^ : V — > Vh- We put 
"projection" in quotes because this need not be the actual orthogonal projection 

with respect to the inner product; indeed, that projection is not generally a 
morphism of Hilbert complexes, since it may not commute with the differentials. 
However, the map ir^ is ^-bounded, surjective, and idempotent. It follows, then, 
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that although it does not satisfy the optimality property of the orthogonal projection, 
it does still satisfy a quasi- optimality property, since 

\\u~ir h u\\ v = inf || (7 - ir h ) (u - v)\\ v < \\I - Tr h \\ inf ||u-w|| v , 

veVh veVh 

where the first step follows from the idempotence of 7T/j, i.e., (I — ir^v — for 
all v G Vh- With this framework in place, the following error estimate can be 
established. 

Theorem 2.10 (Arnold, Falk, and Winther @], Theorem 3.9). Let (V h ,d) be a 
family of subcomplexes of the domain complex (V, d) of a closed Hilbert complex, 
parametrized by h and admitting uniformly V -bounded cochain projections, and let 
(a, u,p) € V k ~ 1 x V k x Sj k be the solution of ([!]) and (a h , u h ,p h ) <E V k ~ l x V k x 
the solution of problem Then 

\W - <rh\\ V + \\u - U h \\ v + \\P-Ph\\ 

<C( inf llcr — t\\ v + inf \\u — v\\ v + inf \\p 
rev*- 1 vev* g ev* 

where fj, = /xjj = sup || (/ — 7r*) r|| . 

||r||=l 

Therefore, if Vh is pointwise approximating, in the sense that inf^gv,^ \\u — v\\ — > 
as h — ¥ for every u £ V, then the numerical solution converges to the exact solution. 

2.5. Improved error estimates. Finally, it can be shown that one can establish 
improved estimates in the VF-norm, subject to a "compactness property." The 
Hilbert complex (W, d) is said to have the compactness property if V k n V fc * is a dense 
subset of W k , and if the inclusion I: V k n V£ W k is compact. Furthermore, 
assume that the family of projections 7T/j is uniformly VF-bounded (rather than 
merely ^-bounded) with respect to h. These properties hold for many important 
examples — notably the L 2 -de Rham complex of differential forms — and allows for 
an abstract generalization of duality-based L 2 estimates (i.e., the Aubin-Nitsche 
trick) to the mixed variational problem. 

The compactness of the inclusion implies that K is also compact, so one may 
define the coefficients 

s = s h = Wi 1 ~ ^h) K\\^ wk wk) , M = Mh = U 1 - n h) P^\\c(w k ,W) > 
V = Vh = max 1 1|(7 - n h ) dK\\ c{wk ^ iWk . j+1) , ||(7 - n h ) A*K\\ c{wk+ ^ wk+j ^ , 

each of which vanishes in the limit as h — > 0. Next, let us denote best approximation 
in the VF-norm by 

E(w) = inf \\w-v\\ , w G W k . 

vev£ 

Then the improved estimates are stated in the following theorem. 

Theorem 2.11 (Arnold, Falk, and Winther g], Theorem 3.11). Let (V,d) be 
the domain complex of a closed Hilbert complex (W, d) satisfying the compactness 
property, and let (Vh, d) be a family of subcomplexes parametrized by h and admitting 
uniformly W -bounded cochain projections. Let (a,u,p) € V k ~ 1 x V k x S) k be the 
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solution of (JlJ and (crhi u hiPh) £ V h 1 x x &h ^ e solution of problem Then 
for some constant C independent of h and (a,u,p), we have 

\\d(a-a h )\\<CE(da), 

\W-o- h \\ < C [E(a) + r,E (da)} , 

\\p-ph\\<C[E(p)+fiE(da)], 
\\d(u- u h )\\ <C(E (du) +v[E (da) + E(p)\) , 

||«-«h|| < C(E(u)+ V [E(du) + E(a)} 

+ {if + S) [E (da) + E(p)} + fiE (P<su)) . 

For typical applications to the de Rham complex, V h k consists of piecewise 
polynomials defined on a mesh. In this case, the order of these coefficients is given 
by rj = 0(h), S = O (ft mm ( 2 ' r + 1 )) , and \x = (h r+1 ), where r is the largest degree 
of complete polynomials in Vu (Arnold, Falk, and Winther [4, p. 312]). 

2.6. Variational crimes. More generally, suppose that the discrete complex Vh is 
not necessarily a subcomplex of V, but that we merely have a VF-bounded inclusion 
map ih '■ Vh V, which is a morphism of Hilbert complexes. Furthermore, given 
the 1^-bounded projection map ir^ : V — > Vh , we require that ^\°i\ — \dyh for 
each k (which corresponds to the idempotence of iTh when ih is simply the inclusion 
of a subcomplex Vh C V). When ih is unitary — that is, when the discrete inner 
product satisfies (u} l: Vh) h = (ih u h,'ih v h) for all Uh,Vh £ V k — then this is precisely 
equivalent to considering the subcomplex ihVh C V. However, if ih is not necessarily 
unitary, we have a generalized version of the discrete variational problem ([2]), stated 
as follows: Find (a h ,u h ,Ph) £ V^ 1 x v£ X Sj k satisfying 

(o-h, Th) h - ( u h, dhTh) h = °) Vt/j £ V*' 1 , 

(3) (dhO-h,v h ) h + (dhUh,d h Vh) h + (Ph,v h ) h = (fh,v h ) h , Vv h £ V k , 

(uh,qh) h = 0, Vq h £Sih- 

The additional error in this generalized discretization, relative to the problem on the 
subcomplex ihVh C V, arises from two particular variational crimes: one resulting 
from the failure of ih to be unitary, and another resulting from the difference between 
fh an d i*,J- 

In Hoist and Stern [24j , we analyze this additional error by introducing a modified 
problem on Vh, which is equivalent to the subcomplex problem on ihVh C V . Define 
Jh = i* h ih, so that for any u hl v h £ W h , we have (ihU hl ihVh) = (ilihUh,v h ) h = 
{JhUh,Vh) h - (The norm \\I — Jh\\, therefore, quantifies the failure of ih to be 
unitary.) This defines a modified inner product on W k , leading to a modified Hodge 
decomposition W% = 05^ © Sj h k 3ft L ' w , where 

tin = 3h I ihz -L h&h} , ~bh x ' w = {v e w£ \ i h v J_ i h 3 k h } . 

Then the subcomplex problem is equivalent to the following mixed problem: Find 
Wh^P'h) e V*- 1 x V* x Sft satisfying 

(J h a h , T h ) h - {Jhu'h, d h T h ) h = 0, Vt,j £ Vjf" 1 , 

(4) {J h d h a'h,Vh)h + (Jhdhu'h,d h v h ) h + (J h p' h ,v h ) h = (i* h f,v h ) h , Vu h £ V% , 

(J h u' h ,q' h ) h =0, Vq'ht&h- 
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The additional error, between the generalized problem ^ and the subcomplex 
problem Q, is estimated in the following theorem. 

Theorem 2.12 (Hoist and Stern [21], Theorem 3.9). Suppose that (<Jh,Uh,Ph) € 
V^ -1 x V k x ft* is a solution to ^ and (a' h , u' h ,p' h ) E V^ -1 x V k x fi'fi is a solution 
to Q. TTien 

IK - <\\v h + IK - «fcllv h + IK -KIU < C - i* h f\\ h + \\I- J h \\ 11/11) . 
Using the triangle inequality, together with the previously stated result of Arnold, 



Falk, and Winther (Theorem 2.10) for the subcomplex problem, we immediately get 
the following corollary. 

Corollary 2.13 (Hoist and Stern [23|, Corollary 3.10). If (a, u,p) E V^ 1 xV k xfi k 
a solution to ([T]) and (crh,,Uh,Ph) € V^ 1 x V h k x $j k is a solution to then 



is 



\\o - ihO-h\\ v + \\u - ihUh\\ v + \\p- ihPhW 

<C( inf lie — t|| v + inf \\u — v\\ v + inf \\p — q\\ v +n inf \\Prgu — v\ 

reihV*- 1 vei h v* q£ihV£ vei h v£ 



where /i is defined as in Theorem 2.10 



This raises the question of how to choose fh € Wh such that ft — > i* h f as h — > 0. 
While fh — i* h f would be the ideal choice, of course, it may be difficult to compute 
the inner product on W , and hence to compute the adjoint it . The following result 
shows that, if : W k — > is any bounded linear projection (i.e., satisfying 
H h o i|F = id W k), then choosing fh = Hhf is sufficient to control this term. 

Theorem 2.14 (Hoist and Stern [H], Theorem 3.11). IfU h : W k -> is a family 
of linear projections, bounded uniformly with respect to h, then we have the inequality 

l|n fc /-*l/|| fc <C(||/- A|| H/II+ inf 11/ -^||). 

Thus, if the family of discrete complexes satisfies the "well-approximating" con- 
dition, and if ||J — Jh|| — > as h — > 0, then it follows that the generalized discrete 
solution converges to the continuous solution. 

3. Semilinear mixed problems 

3.1. An alternative approach to the linear problem. In this subsection, we 
introduce a slightly modified approach to the linear problem, which will be more 
useful in the nonlinear analysis to follow. 

Consider the linear operator L = L ® Pfi : Dl —> W k . Given any u E Dl, we 
can orthogonally decompose u = u + p, where p = P^u and u — u — p. Therefore, 

Lu = Lu + PfiU = Lu + p, 

so given some / E W , solving Lu + p = / is equivalent to solving Lu = /. 
Furthermore, if we define the solution operator K = K © P^, it follows that 

Kf = Kf + P s -J = u+p = u, 

so K is in fact the inverse of L. Thus, L and K establish a bijection between Dl 
and W k . Effectively, by adding P% to each of the operators L and K, we have 
managed to remove their kernel fy k . 
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This approach also sheds new light on the well-posedness of the linear problem. If 
u is a solution to Lu = /, then it satisfies the variational problem: Find u £ V flV£ 
such that 

(5) (d*u, d*v) + (du, dv) + (P^u, P,jv) = (/, v), e v k n V£. 

In fact, the left-hand side is precisely the inner product (u, v) vnv „, which is 
equivalent to the usual intersection inner product obtained by adding the inner 
products for V and V* (Arnold, Falk, and Winther [U p. 312]). Hence, by the 
Riesz representation theorem, a unique solution u = K/ exists, and moreover K is 
bounded. In particular, this variational formulation also illustrates that K is the 
adjoint to the bounded inclusion X: V k C\V^ <-» W k , with respect to this (•, -) VnV * 
inner product, and thus K must be bounded as well. 

Remark 3. While the solutions to the two variational problems ([!]) and ^ are 
equivalent, the mixed formulation is still preferable for implementing finite element 
methods, since one may not have efficient finite elements for the space V k n V fc *. 
We emphasize that this alternative approach is introduced primarily to make the 
analysis of semilinear problems more convenient. 

3.2. Semilinear problems and the abstract Hammerstein equation. Given 
some / £ W k , we are interested in the semilinear problem of finding u, such that 

(6) Lu + Fu = f, 

where F : V k — > W k is some nonlinear operator. Extending the argument from 
the linear case, it follows that this operator equation is equivalent to the mixed 
variational problem: Find (a,u,p) £ V k ~ 1 x V k x S) k satisfying 

(a,T)-{u,dr)=0, Vt £V k ~ 1 

(7) (da,v) + (du,dv) + (p,v) + (F(u + p),v) = (f,v), Vv£V k , 

{u,q)=0, Vq£fi k . 

In the special case where F = 0, this simply reduces to the linear problem. 
Using the solution operator K, the equation Q is also equivalent to 

(8) u + KFu = Kf. 

Equations having this general form are called abstract Hammerstein equations, and 
are of particular interest in nonlinear functional analysis (cf. Zeidler [30j). This 
formulation, which notably appeared in the seminal papers of Amann [TJ [5] and 
Browder and Gupta [S], generalizes certain nonlinear integral equations, called 
Hammerstein integral equations. (In the context of integral equations, the operator 
K corresponds to the kernel operator, or Green's operator.) 

3.3. Well-posedness of the semilinear problem. Before we establish the well- 
posedness of the abstract Hammerstein equation pi), it is necessary to define some 
special properties that a nonlinear operator may have. 

Definition 3.1. The operator A: W k — > W k is said to be monotone if, for all 
u, v £ W , it satisfies (Au — Av,u — v) > 0. It is called strictly monotone if 
(Au — Av, u — v) > whenever u ^ v, and strongly monotone if there exists a 
constant c > such that (Au — Av, u — v) > c \\u — v\\ . 
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Definition 3.2. The operator A: W k — > W k is said to be hemicontinuous if the 
real function t 1— > (A(u + tv) , w) is continuous on [0, 1] for all u,v,w £ W k . 

Theorem 3.3. If F is monotone and hemicontinuous, then the semilinear problem 
^ has a unique solution. Moreover, the problem is well-posed: given two functionals 
f and f , the respective solutions u and u' satisfy the Lipschitz continuity estimate 
||u-u'|| VnV .<||K||||/-/'||. 

The existence/uniqueness portion of the proof is an adaptation of a standard 
argument for Hammerstein equations, when the kernel operator is symmetric and 
monotone on some real, separable Hilbert space (cf. Zeidler [301 P- 618]). 

Proof. Let us define the operator A = I + J£F on V k n V fc * , so that the abstract 
Hammerstein equation ([8]) can be written as Au = K/. Since F is hemicontinuous, 
it follows that A is also hemicontinuous. Moreover, A is strongly monotone with 
constant c = 1, since for any u, u' £ V k D V fc *, we have 

(Au - Au', u - uVnv = H u - u 'llvnv + ( K - ^ u ') , u - u ')vnv* 

ii.. __/ii2 



> u-u' 



Wnv* 

2 

Wnv* 



(Fu-Fu',u u'> 



where the last line follows from the monotonicity of F. Therefore, since A is 
hemicontinuous and strongly monotone, the Browder-Minty theorem [5S] implies 
that it has a Lipschitz continuous inverse A~ l with Lipschitz constant c _1 = 1. 
Hence, there exist unique solutions u = A _1 K/ and u' = A~ 1 "Kf . Finally, by the 
fact that A^ 1 is nonexpansive, these solutions satisfy 

||u - u'|| vnv . < ||K/ - Kf'\\ vnv , < ||K|| ||/ - /'|| , 

which completes the proof. □ 

3.4. Solution estimate for the mixed formulation. Now that we have estab- 
lished the well-posedness of the semilinear problem ([6]), we can use the linear solution 
theory, as developed by Arnold, Falk, and Winther [I], to develop a similar estimate 
for the mixed formulation. This requires placing slightly stronger conditions on the 
nonlinear operator F. In particular, we require F to be Lipschitz continuous with 
respect to the T^-norm: that is, there exists a constant C such that 

||Fu-Fu'||<C||u-u'|| v , 



for all u,u' G V . (Later, in Section 4.5 we will see how this condition can be 



relaxed in case F is only locally Lipschitz.) 

Theorem 3.4. If F is monotone and Lipschitz continuous with respect to the 
V -norm, then the mixed semilinear problem |7| has a unique solution (o~,u,p). 
Moreover, the problem is well-posed: given two functionals f and f, the respective 
solutions (a,u,p) and (a',u',p') satisfy the Lipschitz continuity estimate 

\\<T-0- / \\ v + \\u-u'\\ v + \\p-p'\\<C\\f-f\\, 

where the constant C depends only on the Poincare constant cp and on the Lipschitz 
constant of F. 
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Proof. If u is a solution of the semilinear problem Lu + Fu = /, then it is also a 
solution of the linear problem Lu = g, where g — f — Fu. Therefore, (a,u,p) G 
V k ~ 1 x V k x F) k is the unique solution of the mixed linear problem with functional 
g, and hence of the mixed semilinear problem |7| . 

Now, suppose that u' is the solution to Lu' + Fu' = /', and hence to the linear 
problem Lu' = g' = /' — Fu'. Define u = u u' and g = g — g'; subtracting 
the two linear equations Lu = g and Lu' = g', it follows that Lu = g. Therefore, 
(a, u,p) = (er — a' ,u — vf ,p — p') satisfies the mixed linear problem with functional 
g, so by the well-posedness of the mixed linear problem, we have 

INv + INIv + llPlI <c||ff||, 

where c depends only on the Poincare constant cp. Next, the right-hand side can 
be estimated by 

||s||<||/-/'|| + ||Fu-Fu'|| 

<||/-/'|| + C||u-u / || v <||/-/ / ||+C||u-u'|| VnV ,, 

using the Lipschitz property of F. Finally, applying the previously-obtained estimate 
||u-u'|| vnv . < ||K||||/-/'||, we get \\g\\ < C \\f - f% so finally 

\\a-a'\\ v + \\u-u'\\ v + \\p-p'\\ <C\\f-f'\\, 

which completes the proof. □ 

Remark 4. Note that, in the linear case where F = 0, we can take /' = so that 
(a',u',p') — 0. Then, since g = f and g' = f = 0, we simply recover the usual 
linear estimate \\cr\\ v + \\u\\ v + \\p\\ < c ||/||. 



4. Approximation theory and numerical analysis 

4.1. The discrete semilinear problem. To set up the discrete semilinear prob- 
lem, and develop the subsequent convergence results, we begin by assuming the same 
conditions as in the linear case. Namely, suppose that Vh C V is a Hilbert subcom- 
plex, equipped with a bounded cochain projection 7r^ : V — > Vh- Let Kh : W k — > 
be the discrete solution operator for the linear problem, taking P^f <-> Uh- As with 
the continuous problem, we define a new solution operator = Kh © P$j h and 
consider the discrete Hammerstein equation 

u h + K h P h Fu h = K h P,J. 

Note that this is not simply the Galerkin problem for the original Hammerstein 
operator equation ([8]), since is not just a projection of K onto the discrete space; 
in particular, we generally have 9)\ <£_ fr . 

This is precisely the abstract Hammerstein equation on the discrete Hilbert 
complex Vh, in the sense of the previous section. Therefore, there exists a unique 
solution Uh, and the discrete solution operator Phf H> U/j, Phf H> u' h , satisfies the 
Lipschitz condition 

IK - <\\v h nv : < W K h\\ \\P h (/ - /Oil < l|K h || ||/ - /'|| . 
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Equivalently, this gives a solution to the discrete mixed variational problem: Find 

(a h ,u h ,p h ) e V^" 1 x Vfi x ^ satisfying 

(a h ,T)-(u h ,dT)=0, Vre^f- 1 , 
(9) (da h , v) + (du h ,dv) + (p h , v) + (F (u h +p h ) , v) = (/, v) , Vi> € V£, 

(u h ,q)=0, VgeflJ. 
If F is Lipschitz, then we also obtain an estimate for the mixed solution, 

IK - VhWv + IK - "fcllv + \\Ph-p'h\\ <c h \\f- /I • 

Finally, we remark that when is a family of subcomplexes parametrized by h, 
and the projections 717, '■ V — > Vh are bounded uniformly with respect to h, then the 
constants in these estimates may also be bounded independently of h. 

4.2. Convergence of the discrete solution. We now estimate the error in ap- 
proximating the solution of the mixed semilinear problem ([7| by that for the 
discrete problem Q. Despite the introduction of nonlinearity, we obtain the same 
quasi-optimal estimate as in |Thcorcm 2.10] for the linear problem. 

Theorem 4.1. Let (Vh,d) be a family of subcomplexes of the domain complex (V, d) 
of a closed Hilbert complex, parametrized by h and admitting uniformly V -bounded 
cochain projections, and let (a~,u,p) £ V k ^ 1 x V k x Sj k be the solution of |7| and 
(ch,Uh,ph) € V^" 1 x V£ x fjk £/j e solution of problem ([9]). Then, assuming the 
operator F is Lipschitz with respect to the V -norm, we have the estimate 

Ik - o-fcllv + W u ~ + lb _ Phll 

< C( inf ||cr — r|| v + inf \\u — v\\ v + inf \\p - q\\ v + fi inf \\P<su — v\\ v ), 
Vev* -1 ueVjf gev* t>ey h fc 

where [i is defined as in \Theorem 2.1(\ and where the constant C depends only on 
the Poincare constant cp and the Lipschitz constant of F. 

Proof. Recall that, since [a, u,p) solves the semilinear problem for the functional 
/, it also solves the linear problem for the functional g = f — F(u + p). Let 
(a' h ,u' h ,p' h ) G V^ 1 x Vfi x Sj k be the solution to the corresponding discrete linear 
problem for g. By |Theorem 2.10} this satisfies the error estimate 

Ik - °illy + \\ u - u 'h\\v + \\P-Ph\\ 

<C( inf \\(T — t\\ v + inf \\u — v\\ v + inf \\p — q\\ v + fi inf \\P<su — v\\ v ). 

Vev^ 1 vev* qev£ vev* 

Next, observe that (o~' h , u' h ,p' h ) is also a solution of the discrete semilinear problem 
with functional f = f — F (u + p) + F (u' h + p' h ), since we can just add F (u' h + p' h ) 
to both sides of the equation. However, since the discrete solution operator is 
Lipschitz, we have 

K - a' h \\ v + |K - <\\ v + \\Ph~ P'J <c\\f- /'II 

= C\\F(u+p)-F(u' h +p' h )\\. 

Furthermore, since F is also Lipschitz, 

\\F(u + p)-F(u' h +p' h )\\<C(\\u-u' h \\ v + \\p-p' h \\), 
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which implies 

Ik* - °hlly + W u >> ~ u 'h\\v + \\Ph -p'h\\ < C (|k - a' h \\ v + \\u - u' h \\ v + \\p-p' h \\) ■ 
An application of the triangle inequality completes the proof. □ 

As in the linear case, this implies that if Vh is pointwise approximating in V 
as h — > 0, then (<Jh,Uh,Ph) — > (a,u,p). Moreover, the rate of convergence for this 
semilinear problem is the same as that for the linear problem. 

4.3. Improved estimates. We now establish improved estimates for the semilinear 
problem, subject to the compactness property introduced in |Section 2.5| 

Theorem 4.2. Let (V^d) be the domain complex of a closed Hilbert complex (W, d) 
satisfying the compactness property, and let (V/j,d) be a family of subcomplexes 
parametrized by h and admitting uniformly W -bounded cochain projections. Let 
(cr, u,p) G V k ~ l x V k x fi k be the solution of Q and (a h , u h ,p h ) G V h fc_1 x V% x S) k h 
the solution of problem and assume that the operator F is Lipschitz. Then for 
some constant C independent of h and (a,u,p), we have 

\\d(o- - a h )\\ <C[E (da) + E{u) + E (du) + E[p) 

+ r)E(a) + fiE (P m u)] 

Ik - <rh\\ < C[E(a) + E{u) + E (du) + E[p) 

+ (n + 5 + /i) E (da) + nE {P<su)] 

\\u - u h \\ v + ||p -ph|| < C(E{u) + E (du) + E(p) 

+ n [E(a) + E (da)} + (S + fi) E (da) + pE (P^u)) . 

Proof. As before, let (a' h , u' h ,p' h ) G V^ 1 x V k x ft* be the solution to the discrete 



linear problem with right-hand side functional g = f — F(u + p). Then Theorem 2.11 
gives the improved estimates 

\\d(a-a' h )\\<CE(da), 

\\a-a / h \\<C[E(a)+r 1 E(da)}, 

\\p-p l h \\<C[E(p)+ l iE(da)], 
\\d(u- <)|| <C(E (du) + r, [E (da) + E(p)]) , 

||«-«fc|| <C(E(u) + n[E(du) + E(a)} 

+ (n 2 + S) [E (da) + E(p)\ + pE (P B u)) . 

However, in the proof of Theorem 4.1| we saw that each of the terms ||d (a^ — a' h )\\, 
\\a h -a' h \\, and || u h - u' h \\ v + \\p h - p' h \\ is controlled by 

\\°~h - 0-' h \\y + \\Uh - U' h \\y + \\p h -p'J 

<C(\\u-u' h \\ v + \\p-p' h \\) 
< C(E(u) + E(du)+E(p) 

+ rj [E(a) + E (da)} + (S + fi) E (da) + pE (P<su)) . 

Applying the triangle inequality and eliminating higher-order terms, the result 
follows immediately. □ 
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4.4. Semilinear variational crimes. As first discussed in |Section 2.6[ suppose 
now that Vh is not necessarily a subcomplex of V, and let ih : Vh <^-> V and 
7r;j : V — ► Vh be the W-bounded inclusion and 1^-bounded projection morphisms, 
respectively, satisfying 717, ° ih = idy h ■ Given a discrete functional € TV£ and a 
discrete nonlinear operator Fh : —> , we wish to approximate the continuous 
variational problem ^ by the discrete problem: Find (ah, Uh,Ph) G ^/T _1 x ^/f x &h 
satisfying 

{&h,Th) h - {uhAhTh) h = 0, Vt/j G V^ -1 , 

(dhVh,Vh) h + (dhUh,dh,Vh)h + (Ph,v h ) h 

+ (F h (iih+Ph) , Vh) h = (h,Vh) h ) Vw /i G V^, 
{uh,Qh) h = 0, Mq h &tih- 
For the following error estimate, we define the projection map Py h ■ V — > Vh so that 
ihPv h v is the V^-orthogonal projection of u onto the subcomplex ihVh C V. 

Theorem 4.3. Let(a,u,p) G V^'^x V fc x,f) fe be the solution to ([7| and (a h ,u h ,p h ) € 
Vt~ x V^* x i^^ 6e i/ie solution to (10 1, /jf -FX is Lipschitz, and its constant is 
uniformly bounded in h, then 

\W - ih<7h\\v + ll u - + Hp - 

< C( inf ||cr — r|| y + inf ||u — v\\ v + inf — inf \\P<bu - v\\ v 

+ \\i* h (f-F (u + p)) - (f h - F h P Vh (u + p))\\ h + \\I- J h \\ ||/ — F (u+p)\\), 
where fi is defined as in \Theorem 2.1(\ 

Proof. Suppose (a' h ,u' h ,p' h ) G V^ 1 x V£ x fi^ is the solution to the discrete linear 
problem with right-hand side functional i* h g = i* h (f — F (u + p))- Then, applying 



Corollary 2.13 we have 



Ik - ihO-' h \\ v + \\u - ihu'h\\ v + \\p-ihPh\\ 
<C( inf ||o--r|| v + inf \\u-v\\ v + inf \\p-q\\ v +fi inf \\P<su - v\\ v 

reihV*- 1 v£i h V* qei h V* v<Ei h V£ 

+ \\I-Jh\\\\f-F(u+p)\\). 

Next, observe that (a' h , u' hJ p' h ) also solves the discrete semilinear problem with 
right-hand side functional f' h = i* h (f — F (u + p)) + Fh (u' h + p' h ). Therefore, since 
the discrete solution operator is Lipschitz, we obtain 

\Wh-o-' h \\ Vh + |K - u' h \\ Vh + \\p h -p' h \\ h 

< C |K (f-F (u + p)) - (f h - F h « + p' h ))\\ h 
<C\\i* h (f-F (u+p)) - (f h - F h P Vh (u + p))\\ h 
+ \\F h P Vh (u + p)-F h (u' h +p l h )\\ h . 
Applying the Lipschitz property of Fh to the last term of this expression, 
\\F h P Vh (u+p)- F h (u'h+p'h)\\h < C (\\P Vh u - u' h \\ Vh + \\P VhP -p' h \\ Vh ) 

= C (\\P Vh (u - i h u'h)\\ Vh + \\Pv h (P - ihp' h )\\v h ) 
< C (\\u- i h u' h \\ v + \\p-ihPh\\) , 
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which we have already controlled. Hence, an application of the triangle inequality 
completes the proof. □ 

Clearly, the optimal choice for the functional fh and the operator Fh would be 

fh = i*hf, F h = i* h Fi h . 

In this case, we would obtain 

\\e h (f-F (u + p)) - (f h - F h P Vh (u+p))\\ h 

= \\i* h (F(u+p)-Fi h P Vh (u + p))\\ h 
<C\\(I-i h P Vh )(u + p)\\ v 

<C\ inf \\u-v\\ v + inf \\p - q\\ v , 
\vei h v* q ei h v* J 

which already appears elsewhere in the estimate. Hence, this choice of fh and Fh 
allows the term \\i* h (f - F (u + p)) - (f h - F h (P h u + PhP))\\ h to be dropped. 

However, as noted before, it may not be feasible to take fh = i%,f or Fh = ih^^h, 
since it is often difficult to compute the adjoint i* h to the inclusion. Instead, letting 
n^: W k — > Wh be any bounded linear projection, suppose we choose fh = Hhf 
and Fh = TlhFih, effectively approximating i* h by n^. As in the linear case, this 
choice will give us good convergence behavior, contributing an error that is again 
controlled by other terms in the error estimate. 

Theorem 4.4. Given a family of linear projections Hh - W k — » W*, bounded 
uniformly with respect to h, suppose that fh = Hhf an d Fh = TlhFih, where F is 
assumed to be Lipschitz. Then 

\\i* h (f-F (u+p)) - (f h - F h P Vh (u+p))\\ h < C(\\I - J h \\ ||/ — F (u+p)\\ 

+ inf ||(/- F (u + p)) -0||+ inf Ju-v\\ v + inf \\p-q\\ v ). 

Proof. We begin by using the triangle inequality to write 

\\i* h (f-F (u+p)) -H h (f- Fi h P Vh (u + p))\\ h 

<\\(i* h -TL h )(f-F(u + p))\\ h + \\U h (F (u+p)- Fi h P Vh (« + p))\\ h ■ 
For the first term, we can apply [Theorem 2.14| to obtain 

\\(i* h -Hh)(f-F(u+p))\\ h 

<C(\\I-Jh\\\\f-F(u+p)\\+ inf \\(f-F(u+p))-^\\). 

For the remaining term, we have 
Hllfc (F (u+p)- Fi h P Vh (u + p))\\ h <C\\F(u + p)- Fi h P Vh (u + p)\\ 

<C\\(I-i h P Vh )(u+p)\\ v 

<C\ inf ||ti-w|| v + inf ||p-«|| v J, 

which completes the proof. □ 

Hence, we again get convergence of the discrete solution to the continuous solution, 
as long as the discrete complex is well-approximating and \\I — J/J — > as h — > 0. 
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4.5. Remarks on relaxing the Lipschitz assumption. Our a priori estimates 
for the mixed semilinear problem depended, crucially, on the assumption that the 
monotone operator F was not merely hemicontinuous but Lipschitz. In many 
problems of interest, however, F may be only locally Lipschitz: that is, given 
u e V , there exist constants C, M > (possibly depending on u) such that 
H-Fu- Fu'H < C || u — u'||y whenever ||u — u'||y < M. What can we say about 
well-posedness and convergence when the Lipschitz condition is only local rather 
than global? 

Since [Theorem 3~3 requires only the hemicontinuity of F, we still know that the 



semilinear problem has a unique solution, and that it satisfies 

Hu-u'lUv, <||K||||/-/'||. 

For the mixed problem, though, all we can show is that 

||tr - <r'|| v + ||u - it'll v + \\p - p'\\ <C(\\f- /'|| + \\F u - Fu'||) , 

at which point the proof of |Thcorem 3. 4| requires the Lipschitz condition to continue. 
However, if F is locally Lipschitz at u, then we can still proceed to obtain 

Wa-a'Wy + Wu-u'Wy + Wp-p'W <C||/-/'||, 

as long as ||/ — /'|| (and therefore |ju — u'||y) is sufficiently small. The same holds 
true for the well-posedness of the discrete mixed problem on Vu- 

Now, let us observe how this affects the convergence of the discrete problem. In 
the proof of the a priori estimate, |Theorem 4.1| we had 

\\f-.f'\\ = \\F(u+p)-F(u' h +p' h )\\, 
where (o~' h ,u' h ,p' h ) is the solution to the discrete linear problem with right-hand side 



functional g = f — F (u + p). If Vh is well-approximating in V, then Theorem 2.10 
imples that, by taking h sufficiently small, we can get ||/ — /'|| to be as small as we 
want. Therefore, the error estimates hold as long as h is sufficiently small. 

As an example of how these Lipschitz conditions arise, consider the following 
semilinear elliptic problem on a smooth, connected, open domain f2 C K n : Find 
u G H 1 ^) such that 

(11) -Au + u m = f, 

where m > 1 is an odd integer. Since L — —A is the Hodge-Laplace operator for 
the L 2 -de Rham complex when k = 0, this problem can be expressed within our 
semilinear framework by taking Fu — u m . While F is monotone (since m is odd), 
it does not appear to be globally Lipschitz when m > 1, since the inequality 

(12) \\Fu- Fu'\\ Y <C\\u-u'\\ x , Vfi.u'el, 

cannot be shown to hold for any reasonable choice of the spaces X and Y. 

However, for semilinear scalar problems where both continuous and discrete 
maximum principles are available, it is possible to establish a priori L°° estimates 
on the continuous and discrete solutions. These estimates ensure that the solutions 
both lie in an order interval [u_, u+] n H 1 ^) within the solution space. In other 
words, if u and Uh are the continuous and discrete solutions of the semilinear 



problem ( 11 ), then they satisfy 

U_ < U,Uh < u + . 
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This pointwise control makes it possible to establish (12 1 in this order interval, 
where X — H 1 (O) and Y — L 2 (il). This is precisely the Lipschitz condition that we 
need to apply the framework developed in this paper. In fact, even exponential-type 



nonlinearities can be shown to satisfy the condition ( 12 ) at the continous and discrete 
solutions; see, for example, [TU]- For a discussion of these and related techniques for 
semilinear problems, see |29) . 



While pointwise control of the continuous solution to (11) is always available, 
due to the maximum principle property of the Laplacian, pointwise control of the 
discrete solution is in fact a much more delicate property. Typically, this requires 
placing restrictive angle conditions on the mesh underlying the finite element space. 
In two spatial dimensions, the angle conditions necessary to preserve the maximum 
principle property are achievable with careful mesh generation, even when local 
mesh refinement algorithms in are use. However, in three spatial dimensions, it is 
very difficult to satisfy the required angle conditions, even on quasi-uniform meshes. 

Nevertheless, in the case of sub-critical and critical-type polynomial nonlinearities, 
it is possible to establish a local type of Lipschitz condition by relying only on 
pointwise control of the continuous solution, without requiring pointwise control 
of the discrete solution, and thus avoiding the need for mesh conditions altogether. 
For this class of nonlinearities, one can obtain the following local Lipschitz result. 

Theorem 4.5. Let SI C K.™ for n > 2, and assume that ||w||z,°°(r2) < oo. Let 
F: H l (ST) — > _ff _1 (r2) be a polynomial in u with measurable coefficients defined 
on ft, and whose polynomial degree m satisfies 1 < m < oo for n = 2 and 1 < 
m < fn — (n + 2)/(n — 2) for n > 2. Assume also that u,u' £ i? 1 (f2), and that 
\\u — it'|| ^fi^) < M for some finite constant M . Then 

\\Fu- Fu'\\ H -i {n) < C\\u-u'\\ A ^ n) , 

where C = C(f2, F, \\u\\r,^^),n, m, M). 

Proof. See [5]. □ 



We note that the result in |Theorem 4.5| has a slightly different form than that 
considered above, since F; H 1 ^) — > _ff _1 (f2) rather than — > L 2 (fl). In the 
language of Hilbert complexes, that is, the codomain is given by the dual to V k 
instead of W k . However, as remarked by Arnold, Falk, and Winther [H p. 305], 
the estimates of finite element exterior calculus also apply when the data is given 
weakly as / <E (V k ) , equipped with the sup-norm, and the analysis does not 
change substantially from the / £ W k case (although the solution can no longer be 
interpreted as giving the Hodge decomposition of / in a strong sense). Likewise, 
the results presented here for the semilinear problem also extend to the case of 
weakly-specified data, since the tools of monotone operator theory and abstract 
Hammerstein equations carry over without any significant modification (other than 
the appearance of the sup-norm in place of the M^-norm, where appropriate) . 

Finally, many important problems contain nonlinearities satisfying the assump- 
tions needed to establish continuous and discrete pointwise control, either by sat- 
isfying mesh conditions or by |Theorem 4.5] In particular, these examples include 
the Yamabe problem arising in geometric analysis, and the Hamiltonian constraint 
equation in general relativity. For the three-dimensional case, the leading nonlinear 
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terms for both of these problems have the form 

Fu = au J + bu, 

where a, b € L°°(0). Since m = 5 equals the critical exponent fn = (n 
when n = 3, the nonlinearity satisfies the hypotheses of Theorem 4.5 



f 2)/(n - 
See (231 



-2) 
for 



the derivation of pointwise bounds for both problems, using maximum principles. 



5. Conclusion 

In this article, we have extended the abstract Hilbert complex framework of 
Arnold, Falk, and Winther ^ , as well as our previous analysis of variational crimes 
from Hoist and Stern |24) . to a class of semilinear mixed variational problems. Our 
approach used an equivalent formulation of these problems as abstract Hammerstein 
equations, enabling us to apply the tools of nonlinear functional analysis and mono- 
tone operator theory, and to obtain well-posedness results for both continuous and 
discrete semilinear problems. Additional continuity assumptions on the nonlinearity 
yielded a stronger well-posedness result for mixed problems, as well as a priori 
error bounds for the discrete solution. Despite the addition of nonlinear terms, this 
result agrees with the quasi-optimal estimate of Arnold, Falk, and Winther [3] for 
the linear case, and similarly allows for improved estimates to be obtained under 
additional compactness and continuity assumptions. Likewise, in extending the 
variational crimes analysis in |24j to semilinear problems, we obtain convergence 
results agreeing with the linear case. These last results can also be used to extend 
the a priori estimates for Galerkin solutions to the Laplace-Beltrami equation on 
approximate 2- and 3-hypersurfaces, due to Dziuk [17] and Demlow |15j . to the 
larger class of semilinear problems involving the Hodge Laplacian on hypersurfaces 
of arbitrary dimension. 

At the conclusion of Hoist and Stern [24], several open problems are mentioned, 
including the extension of the Hilbert complex framework to more general Banach 
complexes. While the Hilbert complex framework was again sufficient for the 
analysis of semilinear problems presented here, Banach spaces become necessary 
when dealing with more general nonlinear problems. Banach complexes appear to 
lack much of the crucial structure of Hilbert complexes, particularly the Hodge 
decomposition, whose orthogonality depends fundamentally on the presence of an 
inner product. However, if there is additional structure present in a Banach complex, 
such as a Gelfand-like triple structure (e.g., W C H C W* , where H is a Hilbert 
complex), then it may be possible to generalize the approach taken here. 
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